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Abstract 


In this article, a numerical approximation to the solution of the Newell- 
Whitehead equation (NWE) and Cauchy problem of ill-posed non-linear dif- 
fusion equation have been studied. The presented scheme is obtained by 
using the derivative of the cubic B-spline quasi-interpolation (BSQI) to ap- 
proximate the spatial derivative of the dependent variable and first order 
forward difference to approximate the time derivative of the dependent vari- 
able. Some numerical experiments are provided to illustrate the method. 
The results of numerical experiments are compared with analytical solutions. 
The main advantage of the scheme is that the algorithm is very simple and 
very easy to implement. 


Keywords: B-spline quasi-interpolation; convection-diffusion equation; dif- 
ference schemes. 


1 Introduction 


The use of spline function and its approximation plays an important role for 
the formation of stable numerical methods. Usually, a spline is a piecewise 
polynomial function defined in region, such that there exists a decomposition 
of D into subregions in each of which the function is a polynomial of some 
degree d. Also, the function, as a rule, is continuous in D, together with 
its derivatives of order up to (d— 1). As the piecewise polynomial, spline, 
especially B-spline, have become a fundamental tool for numerical methods to 
get the solution of the differential equations [9, 13, 15, 16, 26]. The numerical 
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solutions of partial differential equations by B-spline quasi-interpolation are 
introduced in [2, 5, 17, 20, 25]. 

Nonlinear equations play an important role in various filed of sciences. 
The world around us is nonlinear, so these kinds of equations arise natu- 
rally in a variety of models from theoretical physics, chemistry, and biology. 
The diffusion equation, one of these nonlinear equations, describes density 
dynamics in a material undergoing diffusion. It is also used to describe pro- 
cesses exhibiting diffusive-like behaviour, for instance the diffusion of alleles 
in a population in population genetics. It has also a great deal of application 
in different branches of sciences which have found a considerable amount of 
interest in recent years [1, 3, 4, 11, 14, 18, 23, 24]. 

Consider the nonlinear Cauchy diffusion equation as the following 


Au = ¢(a,t), «x € (a,b), t>0 (1) 


with initial condition 


u(x,0) = f(x), «x € {a,d (2) 


and boundary conditions of the form 


u(a,t) = aol), u(b.t) = oil), #20 (3) 
Ault) = 32 = 2 ((n(ua,t) tal) 5) (a 


such that «(t)u(x,t) + w(t) is positive [3, 11, 14, 23], a,b are constants, 
go(t), gi (t), K(t), w(t), f(x) and (x,t) are known functions and $(xz,t) be a 
smooth function. 

The Newell-Whitehead equation models the interaction of the effect of the 
diffusion term with the nonlinear effect of the reaction term. For instance an 
equation to describe nearly 1D traveling-wave patterns is put forward in the 
form of a dispersive generalization of the Newell-Whitehead equation. The 
Newell-Whitehead equation is written as: 


Up = Ure taut Bv", x e€ [a,b], t>0 (5) 


where a, are arbitrary constants, n is a positive integer and subscripts x 
and t denote differentiation. 
Initial and boundary conditions are 


u(x, 0) = f(x), TE [a, b] (6) 


v(a,t) = galt), v(6,t) = ga(t), #20 (7) 


where f;(x), g2(t), g3(t) are known functions. The rest of this paper is orga- 
nized as follows. In Section 2, we obtain the numerical schemes using cubic 
B-spline interpolation to solve the nonlinear Cauchy diffusion equation and 
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Newell-Whitehead equation. Some numerical examples are solved to assess 
the accuracy of the technique and the maximum absolute errors will be pre- 
sented in Section 3.The conclusion appears in Section 4. 


2 B-spline quasi-interpolant applied to the Cauchy 
problem and Newell-Whitehead equation 


Assume that an interval J = [a,b] is given, denoted by S4a(X,,) the space 
of splines of degree d and class C¢~! on the uniform partition X, = 
{z; =a+ih,i=0,1,...,n} with meshlength h = (b—a)/n. Let a basis of 
Sa(Xn) be {Byar, j =1,2,...,.n +d} where B; a, is the jth B-spline of de- 
gree d for the knot sequence r := Cy) sane where r_qg=?T_qii =... =7_1 =4@, 
ln =Tnt1 = + =Tnta = band r; =x; 0<i<n. Since the cubic spline 
has become the most commonly used spline and we need the second order 
derivatives we use cubic B-spline quasi-interpolation in this paper. 


From nonlinear differential equation (1) we have 
uz = O(a, t) + K(t) (uz, + UUge) + w(t)Uee (8) 


and from discretizing this equation in time, we get 


2 
ul = (o(eiste) +n(t) (((us}f) tub (uon}f) Holle) (used) Gy 
+ uf 
where u*, (te)* ; (ten)* are the approximation of the values u(z, t), ux(z,t), 
Ura (x,t) at (ai,tk), te = kr, and 7 is the time step. For fixed k, we can get 
the cubic quasi-interpolation as follows [19]: 


n+3 


Qsuk = S° yj (u*) By 3r(a) (10) 


j= 


where u* = u(x,t.) and the coefficient functionals are respectively: 
pii(u") = up, 
p2(u®) = 3% (7u§ + 18ut — Suk + 2us) 


pj(u*) = § (—uf_s + 8uf_. — uf ,) 3<j<nt+l1 (11) 


bn42(u®) = & (Quk_3 —9uk_, + 18uk_,+7uk), 
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Using the de Boor-Cox formula [12, 21], the cubic B-spline basis B;3-(x), 
and his derivatives can be computed. 


For u* € C4(I) we have the error estimate [19] as 


|ju* — Qau*]|,, = O (h*) (12) 


For approximate (ue), (Uen)* by derivatives of the cubic B- pane quasi- 
interpolant (10) up to the order h® we can evaluate the value of u* at 2; 
by: 


nt+3 nt+3 
(Qsuy) =e (Qsur) = 2 wite’ wi). (18) 
We set 
Ur (ug, ut, ’ we us (( i aly. ,(uk)), 
(14) 
UE. = (Cag) 3 (ut) sexe Gtik) 
where 
(uy) = (aul). a) = (suk) 4S 0,1. (15) 
By (15) we obtain 
ut =i pu, vt = Lp, 16 
2 h i ; ce h2 2 ( ) 


where D,, Dz € R(+)*("+1) are obtain as follows: 


“11/6 2 39 18 0 0 
3 12 1 =e 0 0 

1: 2) © Be Sie acs 
O° e378 SH 1s 


ooo°o 
oooo 


Di = : : : : : ; 
1/12 —2/3 0 2/3 -1/12 0 
0 1/12) —2/3 0 2/3 —-1/12 
0 0 1/6 -1 1/2 1/3 
0 0 —1/3 3/2 -3 11/6 
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2 «=f 4 -t oO i 0 0 
fc <a ft & «¢ 0 0 
-1/6 5/3 -3 5/3 -1/6 0 0 0 
0 0 


O0- 1/6 5/2 3 5/8. 176.2 

D2 = ; ; ; 

je B/S. 8-5 1/6 0 
0 1/6. 3/3 $/s -176 
0 0 0 1 -2 1 
0 0 -l 4 -5 2 


0 
0 
0 
0 


From the initial conditions (2) and boundary conditions (3), we can com- 
pute the numerical solution of (1) step by step using the scheme (9) and 
formulas (16). For implementation of this method from (2) we have U® = 
(f (xo), f (21), f(an))? and from (16), (9) and (3) the following algorithm 
is obtained 


oe — (f (xo), f(#1), ug f ORI) 
for k = 0,1,...,m do 

U,* — 4 D,U*; 

Usa — yeDaU*; 

ugt! <= go(te+1); 

for i= 1,2,...,..—1do 


uktl e 7 (Feist) + k(tk) (((u-"),) ty (vs"),)) 
cence (Uae). +ub; 


end 

k 

Ts — G1 (trai); 

k+l k+l ~k+1 ) k+l k+1 , k+1). 
UEr + (ug ,Uz 5 U5 pete are 


Considering a maximum time like T that 0 < t < T we have m = T/T. 


Similarly from discretizing the Newell-Whitehead equation (5), we get 
et ax (Cs) + avk + B(ut)”) +u* (17) 


where v*,(Un2)} are the approximation of the values v(a,t),UVrx(x,t) at 
(xi,tx),t, = kr, and 7 is the time step. For approximation of (vz); , 
in relations (10), (11) and (13)-(16) we set v* = vu(z,t,) and replacing 
v®,V,Ve.,i = 0,1,...,n respectively. Then from the initial conditions (6) 


and boundary conditions (7), we can compute the numerical solution of (5) 
step by step. 
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3 Numerical examples 


In this section, two examples of the nonlinear Cauchy diffusion equation 
and Newell-Whitehead equation are considered and will be solved by B-spline 
quasi-interpolation method. To show the accuracy of the present method for 
our examples in comparison with the exact solutions, the amounts of errors 
is given in some mesh points and we report error norm which is defined by 


n—-1 |, exact _ ee 
i 


1 i 
l= > [u — (18) 


i=l 


For the computational work we select the following examples from [7, 8, 10, 
22). 


Example 1. Let us consider the following nonlinear differential equation 


My (jeu (t +5) = x) 2 - 21-9, (x,t) € [0,1] [0,1] (19) 


which has the exact solution u(x,t) = 2e' +t. In (19) ¢(a,t) = —ft—- 


9,«(t) = ge~*, w(t) = (t+ 5)e*. In Tablel, relative errors at different time 
levels are compared with the relative errors obtained by Zakeri et al. in [10]. 


In Figures 1 and 2 exact and numerical solutions are depicted. 


Example 2. Relative errors at different time levels are compared with 
the relative errors obtained by Nourazar et al. [8]. for Eq. (5) with 
a = 3,8 = —4,n = 3,a = 0,b = 1 and t = 1 in Table 2. The exact so- 


Jee 
3 £ = >: The graph of the 
evbx ye 2 2-¥*) 


exact and numerical solution, are shown in Figures 3 and 4. 


lution of this example is u(z,t) = 


Table 1: Comparison of relative errors obtained from proposed method and method in 
[10]. 


x Relative errors of proposed method Relative errors obtained in [10] 

t = 0.25 t = 0.50 t = 0.75 t = 0.25 t = 0.50 t = 0.75 
0.2 8.0999e-08 6.6631e-08 6.9815e-08 2.7500e-08  4.7700e-09  3.4500e-08 
0.4 1.0121e-07 9.3251e-08  1.0112e-07 4.8100e-07  4.3300e-07 — 2.8500e-07 
0.6 8.1254e-08 8.1761e-08 9.1375e-08 2.2800e-06 2.2700e-06 1.9100e-06 
0.8  4.4684e-08  4.7547e-08 5.4249e-08 6.7100e-06 6.7700e-06 —5.8700e-06 


jel;  6.4079e-08 5.9977e-08 6.5621¢-08 = = = 


From the test examples, we can say that the BSQI scheme is feasible and 
the accuracy is better than the multi-quadric quasi-interpolation (MQQI) 
method [6]. Moreover, MQQI method has very close relation to the shape 
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Table 2: Comparison of errors of Example 2 with the errors obtained in [8].(h = 0.02,7 = 
0.0001) 


x Relative errors of proposed method Relative errors obtained in [8] 
#=0.1 t=0.15 t=0.2 t=1 t=0.1 t=0.15 #=0.2 

0.2 4.7533e-06 6.0414e-06 6.6041e-06 5.2440e-07 4.9987e-06  5.6384e-05 3.1193e-04 

0.4  6.8110e-06 8.4592e-06 9.1422e-06 7.1097e-07 6.3997e-06 6.8760e-05 3.6460e-04 

0.8  4.7680e-06 5.4354e-06 5.6195¢-06 4.0217e-07 3.6819e-06 3.7324e-05 —1.8700e-04 


lel,  4.8486e-06 5.8747e-06 6.2722e-06 4.7741e-07 = - = 


Figure 1: The exact solution of Example 1 for h = 0.02,7 = 0.00001 


parameter c in MQ. In fact, the choice of the shape parameter is still a pen- 
dent question. Furthermore, the MQQI is required to calculate derivatives 
of MQ quasi interpolant once for all, which is not easy to compute when h is 
small. Although the accuracy of BSQI is not better than that of other meth- 
ods, we know that, at each time step, the complexity of BSQI is lower than 
theirs. The proposed method is an acceptable and valid scheme. Moreover, 
it can be implemented very easily. 


4 Conclusions 


In this article, we have applied the cubic B-spline quasi-interpolation 
method for solving the nonlinear Cauchy diffusion problem and Newell- 
Whitehead equation. The results have been compared with the exact solu- 
tions and demonstrated the good performance of the schemes. This method 
offers several advantages in reducing computational costs. On the other hand, 
this method is very simple to apply and to make an algorithm. Thus, this 
method may be reckoned as a simple and accurate solver for PDEs and it is 
worthy to note that this method can be utilized as an accurate algorithm to 
solve linear and nonlinear functional equations arising in physics and other 
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Figure 2: The numerical results of Example 1 for h = 0.02,7 = 0.00001 


Figure 3: The exact solution of Example 2 for h = 0.02,7 = 0.0001 


fields of applied mathematics. The computations associated with the exam- 
ples in this article were performed using MATLAB R20138a. 
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